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ABSTRACT 

We discuss recent achievements in light scalar mesons spectroscopy through 
amplitude analysis of charm particle decay and its consequences. The high 
statistics clean samples of charmed mesons, in addition to it's definite J p and 
mass, is turning these decays into a new important environment to study light 
meson physics. We give special attention to the scalar sector favored by a high 
coupling to charm. 



1 Introduction 

Most mesons are well understood in the context of the quark model. This is not 
true just for the scalar sector, that for long has been a source of controversy. 
For the latest 30 years a great deal of experimental effort have been made but in 
many cases the experimental results do not converge to compatible outputs. In 



production experiments, observing the scalar resonant states is difficult due to a 
large contribution of the non-resonant background. Moreover, light scalars are 
too numerous within a relatively short mass interval. The absence of a unique 
analysis procedure, especially regarding the interference of the background with 
the resonances results in conflicting measurements. In this sense the wider and 
lighter the state the worse, and this is the case of the two states c(500) and K 
that will be discussed here. 

The theoretical interpretation of the lightest scalar meson have also been 
unclear. Many objects - glueballs, KK molecules, multiquark compact states 
- are expected to populate the area. In a recent review article F. Close and 
N. A. Tornqvist 0) discuss the scalars from both experimental end theoretical 
points of view. They suggest that the scalars be organized not in one but 
two nonets. The "standard" quark model qq nonet, distorted by a glueball 
predicted by lattice QCD, is enough to explain the region above ~ I GeV. This 
nonet is composed by a (~ 1400), / (1370),Jf (1430), / (1500) and /o(1710). 
The states below - lGeV - / (980), a (980) and possibly ct(500) and k - by 
arguments based on QCD attractive forces in S-wave, would also form a nonet. 
The interpretation of such nonet would be more complex, of the type meson- 
meson. In the article, they stress the importance of charm decay as "opening 
up a new experimental window for understanding light meson spectroscopy and 
specially the controversial scalar meson which are copiously produced in these 
decays" . 

The use of charm decay is an alternative to the traditional production 
experiments to study lighter resonant particles. It was made possible by the 
large clean samples of charm now available and which are attributed to mankind 
effort but also, and more important, by what can be seen as nature's gifts: 

1) Non-leptonic charm decays are preferentially two body or quasi- two 
body, through the formation of intermediate resonant states. With a small 
non-resonant (NR) component one avoids having to deal with model dependent 
amplitudes and its interference with the resonant amplitudes. This feature is 
particularly important for the wide scalars because their amplitudes and the 
NR (usually taken as constant) can become very similar. We return to this 
issue later in the text. 

2) Charm couple strongly with scalars. This empirical fact is present 
in the decays discussed here where the contribution to scalars dominates all 



processes. 



3) We are dealing with a well denned initial state; (charm meson mass, 
M and J p ) 

The first light resonant parameters extracted from charm data use D me- 
son decays to three charged pseudo-scalars. In the following session we sum- 
marize the Dalitz plot amplitude formalism used. The parameterization of the 
overall amplitude consists of a coherent sum of all contributions (NR flat and 
resonant as relativistic Breit-Wigners modulated by angular momentum con- 
servation functions and form factors) with complex coefficients obtained from 
the fit. The magnitude of the coefficients are proportional to each relative 
contribution and the phases accommodate in an effective way the final state 
interactions. It should be pointed out as the a great pro of these analysis the 
very good description that such simple model give to the data. Alternative 
descriptions for the decay amplitude are being tried, for example Focus exper- 
iment is using K-matrix to parameterize the light scalars contribution, but we 
shall not discuss those, still preliminary, results. Also expected for the near 
future are Focus four-body amplitude analysis. 

2 Three body hadronic decay formalism 

We describe here the analysis procedure used by the E791 collaboration in 
their D + and D+ — > 7r + 7r~7r+ IP i) and D + — > K~7r + n + B), from which was 
measured masses and widths of <j(500), /o(980), /o(1370), k and #£(1430). 

The decay of a scalar hadron of mass M into 3 spinless daughter particles 
is completely specified with two degrees os freedom, conveniently chosen as 
two Dalitz plot variables, m\ 2 and m| 3 . The Dalitz plot density distribution 
is proportional to the invariant decay amplitude A squared and reflects the 
dynamics of the decay process. A simple analytical model for A is given by a 
coherent of all intermediate states contributing, resonant or not: 



The parameters a give the various relative contributions and the phases 
5 are accounts for final state interactions. The non-resonant amplitude Anr is 
represented by a constant. Which is a reasonable assuption if we imagine that 
it is dominantly 5-wave, in any case, the impact of this choice in the results is 
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reduced due to a small NR contribution observed. The decay through resonant 
intermediate states, j (that then decay to the observed k and I hadrons), are 
viewed as s-channel processes where the resonance plays the role of massive 
propagators, represented by relativistic Breit-Wigner functions, BWj [] . Mo- 
mentum dependent form factors, Fjj and Fr describe the non-pointlikc nature 
of the D meson and the resonance respectively and depend on the resonance 
spin, J and the radii of the relevant mesons. The angular momentum conserva- 
tion is taken care by the function Mj- Each resonant amplitude, Aj is written 
as: 



Aj = BWj xF D xF R x Mj (2) 
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Above mf; is the the Dalitz plot variable, i.e. the invariant mass of the two 
hadrons forming a spin-J resonance. Detailed expression of all the above func- 
tions are found in the references . Combinatorics background, detector 
efficiency and Bose symmetrization are also considered in a maximum likeli- 
hood fit to the data to extract the parameters. In most cases the masses, mo, 
and widths, Tq of the resonances are fixed by values listed in PDG. Only for 
new or poorly measure states they are allowed to float in the fit. 



1 For the /o(980)7r + E791 uses a coupled channel Breit-Wigner, following 
the parameterization used by the WA76 Collaboration , 

BWm = ^-^ + L o(r7r + rvr with 
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Figure 1: Dalitz-plot of the individual resonant contributions to the D + — > 
7T _ 7r + 7r + decay. 

Figure |l| shows , as an example, Monte Carlo simulation Dalitz-plot of 
each individual resonant state that contributes in the D + — > 7r~7r + 7r + decay. 
Notice that being coefficients and the individual amplitudes complex quantities, 
interference effects will take place when all the pieces act together. 

The above described produces a quite stringent model and obtaining ac- 
ceptable fits is usually a difficult task. To access the quality of each fit a fast 
Monte Carlo (MC) program was developed which produces Dalitz plot event 
densities accounting for signal and background PDF's, including detector effi- 
ciency and resolution. Comparing the MC density distribution generated using 
parameters extracted from a given fit with that for the data, it is produced 
a x 2 distribution. When comparing two possible models, the best discrimi- 
nating power test requires ensembles of Monte Carlo "experiments" . In the k 
discussion below we illustrate the technique. 



3 The /qS resonances 



We have chosen to start the results sessions with the /o because of the extremely 
clear signal seen in the _D+ — > n + iT — n + decay, figure ^|. There are dozens of 
measurements listed in PDG 0) for this state and they converge to a reasonably 
well defined mass, 980 ± 10 MeV but the estimation for width is from 40 to 100 
MeV. 

For the best E791 fit 5 resonant channels contributes significantly: p° (770)tt + , 
p°(1450)7r+, / (980)7T+, / 2 (1270)7r+, and / (1370)7r+, plus a NR that con- 
tributed with a fraction of only 0.5±0.2%. The fit have y^jdoj = 71.8/68 with 
a confidence level of 35%. The dominant contributions comes from /o(980)7r + , 
56.5 ±6.4%, and /o(1370)7r + , 32.4±7.9% from which they measure: mj (g 80 ) = 
977 ± 3.6 MeV/c 2 , g^ = 0.09 ± 0.01, g K = 0.02 ± 0.05, m fo{1370) = 1434 ± 20 
McV/c 2 and r /o(1370) = 172 ± 33 MeV/c 2 . There is no evidence of a third 
higher mass scalar state, Jo (1500). 

Previous production experiments claim a large contribution of the KK 
channel by estimation a large value of the parameter gK ^ ® ■ These results 
do not agree with the value measured by E791 that, in fact, finds that a simple 
Breit-Wigner || is sufficient to represent the data. From the small but clear 
preliminary signal BES collaboration measure m/ = 0.980 ± 0.009 GeV and 
T f — 0.045 ±0.30 GeV & . In figure || we compare various results and conclude 
that in charm decay the resonance /o(980) presents a narrower signal. 

4 The cr(500) 

The 1999 the Workshop on Hadron Spectroscopy H) devoted one entire session 
to the meson cr(500), " What do we know about the cr?". By that time E791 
had not published their observation of the meson a in the decay of D + — > 
7r + 7r~7r + E§). Today charm decay is viewed as a tool for studying the light 
spectroscopy. 

In figure ^ we show the projection of the E791 Dalitz plot where a clear 
peak at low mass can be seen. The figure compares the best fit achieved with 
all possible well established resonances available at the time, a) (Fit 1), to 
their solution including a low mass scalar state, b) (Fit 2). They measured 
= 478 ± 29 MeV/c 2 and r CT = 324 ± 46 MeV/c 2 and the inclusion of the 
state took them from an unacceptable solution of X 2 /dof = 254/162 with a 




Figure 2: E791 Df — > ir + ir — ir + Dalitz plot and m + _ projection 



confidence level less than 10" 5 to a very good fit with \ 2 /dof = 138/162 and 
confidence level of 90%. In Fit 1 the NR contribution is dominant with 38% 
and in Fit 2 following the trend of charm decays it dropped to 7.8% whereas 
the er7T + dominates with 46%. E791 perform a series of alternative fits and 
tests to be sure that no other model would as well describe the data. 

Studying the channel J/^ — > u>7r + 7r _ BES experiment observe a signal 
of the mesons a and measure their parameters 0); m a = 490^™ MeV/c 2 and 
T a = 282^5Q MeV/c 2 . In figure || we show their final fit and scans for the a 
mass and width measurements. 



5 The k 



As a last result we consider the k resonance observed in E791 Cabibbo favored 
decay D + — » K~tt + it + . Despite the large statistics available in this channel, no 
previous experiment 0) 111) have been able to provide a convincing explanation 
for this decay. E687 EH) best model for their sample of almost 9000 events have 
X 2 /dof — 87/29. The solution have the NR contribution dominating with 



Figure 3: Comparison of several /o(980) results. 



a fraction of 99%. Large interference pattern is produced with all fractions 
summing 147%. In contrast to the a where a clear bump is seen in the m + _ 
projection, no evidence for a missing piece in the low m 2 K _^ + region can be 
easily observed (figure ^b). On the other hand the very large statistics and 
small number of possible intermediate states provide strong evidence of the 
need for an extra low mass wide resonance contributing for the decay. 

The first approach tried by E791 @) was to include all established states; 
K*{892)n+, i^(1430)7r+, /? 2 *(1430)tt + , if*(1680)7r+ plus a NR (they have stud- 
ied also the possibility of a non-flat NR contribution without success). In gen- 
eral the result agree with previous studies including the bad quality of the 
fit. At this point they measured mass and width of the scalar Xq(1430) to be 
1416±27 and 250 ± 21 MeV/c 2 respectively, which agree with PDG values. 

Next they include an additional scalar named "k" for which they measure 
respectively mass and width of 797± 47 and 410± 97 MeV/c 2 . In the same fit 
the parameters relative to i^Q (1430) were measured to be 1459±9 and 175±17 
MeV/c 2 . This new model describes very well the data with a x 2 /dof = 46/63, 
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Figure 4: El 91 D + — > 7r + 7r — ir + m + _ projections for data (error bars) and 
model (solid line). Shaded area is the background, a) solution for Fit 1, and b) 
Fit 2. 

confidence level of 95%. The contribution of the new state is dominant with 
48± 12% of the total fraction and the NR contribution is of 13±7 %. The sum 
of all fractions dropped from 134% in the fit without k to 88% indicating a 
smaller degree of interferences. 

In the amplitude analysis described here one do not measure directly a 
Breit-Wigner phase, instead it is assumed and consistency tests have to be 
performed to verify if alternative models are able to represent the data. This 
was done for this analysis; a toy-model (consisting of a Breit-Wigner amplitude 
without a phase variation), a vector and a tensor alternative models were tried, 
none producing a satisfying solutions. When comparing two possible models A 
and B the best figure of merit to be used is given by the Neyman-Person EH) 
lema: Awa,b = ~^{lnCA ~ InCs), where £a,b are the likelihood for a given 
set of events calculated with the parameters obtained from fits to the data with 




models A and B. Ensembles of 1000 Monte Carlo "experiments" were generated 
with the parameters extracted from fits to the data for the two models A and B. 
For each of such experiments the quantity Awa,b was calculated and plotted 
in figure f^, where the model with the scalar k is compared with model without 
it; with the toy-model k or with the vectorial k. The discriminating power 
of the exercise is obvious by the separation of the distributions. The value of 
Awa,b for the data is signed by the triangle showing the clear preference of 
the data to the model of the scalar k. 



6 Conclusion 



In conclusion, charm decay is, as said by F. Close and Tornqvist, a new window 
for studying light mesons. We discussed some of its remarkable contributions. 
With the already available data sets and the richness of charm decays, we hope 
to see in the near future confirmations and better measurements for some of 
the resonances, like the k. 
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